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METHOD OF ESTIMATING THE SPATIAL VARIATION OF MAGNETIC 
RESONANCE IMAGING RAD I OFREQUENC Y (RF) SIGNAL INTENSITIES 
WITHIN AN OBJECT FROM THE MEASURED INTENSITIES IN A UNIFORM 
SPIN DENSITY MEDIUM SURROUNDING THE OBJECT 

5 

Field of the Invention 

The present invention relates to a method of estimating the 
10 spatial variation of magnetic resonance imaging 
radiof requency (RF) signal intensities within an object from 
the measured intensities in a uniform spin density medium 
surrounding the object, such as for the subcutaneous fat 
surrounding the liver in an axial image. 

15 

Background of the Invention 

The clinical applications of nuclear magnetic resonance 
imaging (MRI) are progressively expanding into the 

20 quantitative measurement of the physical and chemical 
properties and processes of the human body. Where 
quantitative information is extracted from the image data, 
it is important to account for the variety of instrumental 
effects ttiat may perturb the parameter being measured, as 

25 well as any erroneous effects presented by the object being 
imaged. Of particular importance are the multiplicative 
inhomogeneities contributing to the spatial variation in 
radio frequency (RF) signal intensity within the image. Non- 
uniformities in magnetic resonance signal intensity are 

30 caused by both instrumental and object effects, as perceived 
by J. G. Sled and G. B. Pike, in "Standing-Wave and RF 
Penetration Artifacts Caused by Elliptic Geometry: An 
Electrodynamic Analysis of MRI", IEEE Transactions on 
Medical Imaging, Vol. 17, No. 4, pp. 653-662, 1998 (ref). 

35 Instrumental effects include inhomogeneous RF excitation, 
non-uniformities in receiving coil sensitivity, and gradient 
field eddy currents, whilst measurement object effects 
include variable RF penetration and standing wave effects. 
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These intensity variations are detrimental to the meaningful 
comparison of image intensities in separate parts of the 
image, and must be corrected for if accurate quantitative 
information is to be obtained. Methods that can correct for 
5 the effects of the RF signal intensity gradient are thus an 
important area of development in quantitative MRI. 

A variety of pre- and post-processing techniques have been 
developed with the aim to correct for the spatial variation 

10 in RF signal intensity within the object being measured, as 
discussed by L. Q. Zhou, Y. M. Zhu, C. Bergot, A.-M. Laval- 
Jeantet, V. Bousson, J.-D. Laredo, and M. Laval- Jeantet , in 
"A method of radio-frequency inhomogeneity correction for 
brain tissue segmentation in MRI" , Computerized Medical 

15 Imaging and Graphics, Vol. 25, pp. 379-389, 2001. Pre- 
processing techniques have focussed on the measurement of 
homogeneous phantoms prior to measurement of the object 
under study to estimate the inhomogeneities in the RF field. 
However, these techniques do not typically account for the 

20 inhomogeneities in RF penetration presented by the object 
being imaged. Post-processing techniques attempt to rectify 
this situation by estimating the decay profile of the RF 
signal intensity within the object from analysis of the 
object data set itself. This enables both correction of 

25 object independent effects, such as bias in the RF field, as 
well as object dependent effects, such as signal intensity 
attenuation within the object. The assumption behind most 
post-processing techniques is that the non-uniformity in RF 
signal intensity may be attributed to a low spatial 

3 0 frequency component throughout the image. A number of 
approaches to estimating RF inhomogeneities are thus based 
on intensity correction schemes that employ image smoothing. 
Homomorphic filters are predominantly used in this situation 
as they account for the low frequency signal intensity 

35 components without altering structural boundaries. Other 
post-processing methods are based on structural 
classification techniques and intensity surface 
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interpolation. 

One post-processing technique for the estimation of RF 
signal intensity variations within the human head is worth 
5 particular reference. The method by B. M. Dawant, A. P. 
Zijdenbos, and R.A. Margolin, in "Correction of Intensity 
Variations in MR Images for Computer-Aided Tissue 
Classification", IEEE Transactions on Medical Imaging, Vol. 
12, No. 4, pp. 770-781, 1993, consists of interpolating a 

10 thin-plate spline surface to reference points within the 
white matter of the brain and to use the interpolated 
surface as an estimator of the RF coil profile. The 
availability of reference points on the periphery of the 
images is crucial for obtaining good correction surfaces, so 

15 intensity values of edge points are extrapolated from the 
interior reference points. To reduce the sensitivity of the 
technique to the mislabelling of reference points selected 
by the user, the number of reference points may be increased 
by a tissue classification technique. The thin-plate spline 

20 surface is then fitted to the reference points by the method 
of least-squares rather than by interpolation. 

The main limitation of the intensity surface correction 
technique described above is that there is no measured data 

25 with which to justify the extrapolated intensity reference 
values at the periphery of the head, given that the signal 
intensity gradient is not known at these peripheral points. 
Further, there is no physical basis to the selection of the 
thin plate spline surface to model the RF signal intensity 

30 gradient. As such, the thin-plate spline surface may not 
adequately account for local variations in RF attenuation 
throughout the object being imaged. 

It is to be understood that, if any prior art publication 
35 is referred to herein, such reference does not constitute 
an admission that the publication forms a part of the 
common general knowledge in the art, in Australia or any 
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other country. 

In the claims which follow and in the preceding description 
of the invention, except where the context requires 
5 otherwise due to express language or necessary implication, 
the word "comprise" or variations such as "comprises" or 
"comprising" is used in an inclusive sense, i.e. to specify 
the presence of the stated features but not to preclude the 
presence or addition of further features in various 
10 embodiments of the invention. 

Summary of the Invention 

It is an objective of the present invention to provide a 
15 method of estimating the spatial variation of magnetic 

resonance imaging radiof requency (RF) signal intensities 
within an object from the measured intensities in a uniform 
spin density medium surrounding the object. Once the 
estimation has been obtained it may be used in a variety of 
20 ways to compensate for the spatial variation in RF signal 

intensity within the object or to simulate image intensities 
subject to the RF signal intensity variations. 

According to the present invention there is provided a 
25 method of estimating the spatial variation of magnetic 

resonance imaging radio frequency (RF) signal intensities 
within an object from measured RF signal intensities of a 
uniform spin density medium surrounding the object, said 
method comprising : 
30 acquiring a magnetic resonance image of an object 

bounded by a medium which is of essentially uniform spin 
density, on the length scale of resolution of the image; 



35 



formulating a semi-empirical mathematical model of the 
spatial variation in RF signal intensity within the said 
object; and, 
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fitting said model to selected measured RF signal 
intensities of said surrounding medium to obtain an estimate 
of the spatial variation of magnetic resonance imaging 
radiof requency (RF) signal intensities within said object. 

5 

Preferably said step of formulating said semi-empirical 
mathematical model of the spatial variation in RF signal 
intensity within the object comprises: 

locating a plurality of points in a plane of said image 
10 which are notionally considered to act as apparent receivers 
of RF signals (hereinafter referred to as "RF receiver 
points' 7 ) ; and, 

formulating a spatial profile of RF signal intensity 
15 within said object in said plane of said image relative to 
the said RF receiver points. 

Preferably said RF receiver points are located on said 
surrounding medium. 

20 

Preferably said spatial intensity profile is formulated in a 
manner which provides a concentric reduction in RF signal 
intensity with increasing distance from said RF receiver 
points . 

25 

Preferably said reduction in RF signal intensity is 
formulated as having the same rate of reduction in signal 
intensity with increasing distance from all said RF receiver 
points . 

30 

In one form of said model, said spatial profile is 
formulated as having an exponential reduction in signal 
intensity with increasing distance from said RF receiver 
points . 

35 

However in an alternate embodiment said spatial profile is 
formulated as having a reduction in signal intensity with 
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the reciprocal of the distance from said RF receiver points 
raised to a chosen power. 

Preferably said semi-empirical mathematical model is 
5 formulated as : 

n= 1 

where x and y are image coordinates, I(x,y) is the model of 
the spatial intensity profile of the RF signal, N is the 
number of RF receiver points, I n is the estimated signal 

10 intensity at the n th RF receiver point and is to be 

determined by the fitting procedure, R is the rate of decay 
in signal intensity with distance from said RF receiver 
points and is to be determined by the fitting procedure, x n 
is the x position of the n th RF receiver point, y n is the y 

15 position of the n th RF receiver point, and f is the degree 
of ellipticity to the concentric decay in RF signal 
intensity with increasing distance from said RF receiver 
points and is to be determined by the fitting procedure. 

20 Preferably said selected measured RF signal intensities 

comprise pixels of maximal signal intensity in respective 
radial line segments of an image of the surrounding medium 
(hereinafter referred to as "selected pixels'') . 

25 Preferably said method further comprises connecting said 
selected pixels by a line. 

Preferably said method further comprises deriving from said 
trace local maxima in signal intensity. 

30 

Preferably when said magnetic resonance image is acquired in 
an axial plane said RF receiver points are located at the 
local maxima in signal intensity of the said selected pixels 
from said surrounding medium. 
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Preferably when said magnetic resonance image is acquired in 
said axial plane the said local maxima in signal intensity 
of the selected pixels are no closer than 12 pixels. 

5 Preferably said local maxima in signal intensity determined 
from the selected pixels serve as starting values for the 
signal intensities I n of the said RF receiver points in the 
said fitting procedure. 

10 Preferably said fitting procedure includes deriving an 
initial estimate of the rate R of decay of signal intensity 
with distance from said RF receiver points determined by: 

locating local minima amongst said selected pixels in a 
manner consistent with determination of said local maxima; 

15 

for each RF receiver point locating two nearest local 
minima and fitting an exponential signal intensity decay 
curve as a function of the distance from the said RF 
receiver point to the signal intensities of the two local 
20 minima constrained to the signal intensity at the position 
of the said RF receiver point to determine an estimate of 
the rate R n of decay in signal intensity with distance away 
from the said RF receiver point; 

25 calculating the average of the estimates of the rate of 

decay R n in signal intensity with distance from each RF 
receiver point which is to be used as an initial estimate 
of R. 

3 0 Preferably when the magnetic resonance image is acquired in 
the axial plane the elliptical contours of the decay in the 
RF signal intensity from a given RF receiver point have 
their minor axes approximately normal to a plane of an RF 
receiver coil element of an MRI machine acquiring said image 

35 of said object in closest proximity to the given RF receiver 
point . 
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Preferably when the magnetic resonance image is acquired in 
the axial plane the angular alignment of each of the 
elliptical contours of the decay in the RF signal intensity 
from the RF receiver points is determined relative to said 
5 trace. 



Preferably the angular alignment of the major axes of the 
elliptical contours of the decay in the RF signal intensity 
from the RF receiver points is determined relative to said 
10 trace by: 

centering a window kernel on each RF receiver point and 
finding the two furthest positions from said RF receiver 
point at which said trace of the selected pixels passes 
through the perimeter of the window kernel; and, 

15 

calculating the angle of the line bisecting said two 
furthest positions relative to the coronal plane. 

Preferably said window kernel which is centered on said RF 
20 receiver points to determine the angular alignment of the 
elliptical contours of the decay in the RF signal intensity 
from said RF receiver points is a square window kernel not 
less than 13 pixels in width. 

25 Preferably when the magnetic resonance image is acquired in 
an axial plane said model is formulated as: 

T/ N ^> T - RJ (x-x n ) 2 [l / /. cos 2 d n + /sin 2 6„ 1 (y-y n ) 2 [l / ~ sin 2 0 n +f cos 2 7„ >( x-x n )(y-y n )[l / /.sin 28 n -/sin 2 0„ ] 

where x and y are image coordinates, I(x,y) is the model of 
the spatial intensity profile of the RF signal, N is the 

30 number of RF receiver points, I n is the estimated signal 
intensity at the n th RF receiver point and is to be 
determined by the fitting procedure, R is the rate of decay 
in signal intensity with distance from the RF receiver 
points and is to be determined by the fitting procedure, x n 

35 is the x position of the n th RF receiver point, y n is the y 
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position of the n RF receiver point, f is the degree of 
ellipticity to the concentric reduction in RF signal 
intensity with increasing distance from the RF receiver 
points and is to be determined by the fitting procedure, and 
5 Q n is the angular alignment of the major axis of the 
elliptical contour to the decay in the RF signal intensity 
from the n th RF receiver point with the coronal plane. 

Preferably said model spatial intensity profile I(x,y) is 
10 fitted to the image intensities of the selected pixels from 
the medium bordering the object for optimum parameterisation 
of J n , R and f to obtain an estimate of the spatial 
variation in RF signal intensity throughout the object. 

15 Preferably said method further comprises determining a RF 
spatial attenuation profile with which to rescale the 
magnetic resonance image intensities within the said object 
and thus minimise the spatial variation in RF signal 
intensity throughout the object. 

20 

Preferably said RF spatial attenuation profile is determined 
by dividing the estimated RF signal intensity profile I(x,y) 
by the minimum or maximum value of I(x,y) within the region 
occupied by the object. 

25 

Preferably the spatial variation in RF signal intensity 
throughout the object is be minimised by dividing the image 
intensities within the object by the RF spatial attenuation 
profile . 

30 

Preferably said method further comprises, when the magnetic 
resonance image is a spin echo image, providing an estimate 
of the image intensities within the object at zero echo time 
by dividing the estimated signal intensities within the 
35 object at the given echo time by a percentage remaining 
signal intensity expected within the surrounding medium at 
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the given echo time and then dividing this result by a ratio 
of hydrogen proton spin density expected of the bounding 
medium relative to hydrogen proton spin density expected of 
the object itself. 

5 

Preferably said method also comprises calculating transverse 
relaxation rates within the object from a series of spin 
echo images acquired of the object and surrounding medium 
at different echo times. 

10 

Preferably when the magnetic resonance image is a spin echo 
image, said method further comprises using the estimates of 
the image intensities within the object at zero echo time in 
the calculation of transverse relaxation rates within the 
15 object where the estimates are determined from a spin echo 
time image of shortest duration in the spin echo image 
series . 

Preferably when the magnetic resonance image is a spin echo 
20 image, the transverse relaxation time of said surrounding 
medium is significantly greater than the spin echo time at 
which the image is acquired such that there is negligible 
decay in signal intensity throughout the medium from zero 
echo time. 

25 

Preferably when the object under consideration is the liver 
said surrounding medium is a layer of subcutaneous fat 
surrounding the abdomen. 

30 Preferably when transverse relaxation rates are calculated 
within the liver from a series of spin echo images the 
subcutaneous fat is used to estimate image intensities 
within the liver at zero echo time to assist in the 
calculation of the transverse relaxation rates. 

35 

According to the present invention there is also provided a 
method enhancing a magnetic resonance (MR) image of an 
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object comprising: 

obtaining an initial MR image of the object; 

5 estimating the spatial variation of magnetic resonance 
imaging radio frequency (RF) signal intensities within the 
object from measured RF signal intensities of a uniform 
spin density medium surrounding the object by: 

10 acquiring a magnetic resonance image of an object bounded 
by a medium which is of essentially uniform spin density, 
on the length scale of resolution of the image; 

formulating a semi-empirical mathematical model I(x,y) of 
15 the spatial variation in RF signal intensity within the 
said object; 

fitting said model to selected measured RF signal 
intensities of said surrounding medium to obtain an 

2 0 estimate of the spatial variation of magnetic resonance 

imaging radiof requency (RF) signal intensities within said 
object; and, 

utilising the estimated spatial variation in RF signal 
25 intensity to rescale the initial MRI image by dividing the 
signal intensities of the image object by the RF signal 
intensity profile I(x,y) . 

Brief Description of the Drawings 

30 

An embodiment of the present invention will now be described 
by way of example only with reference to the accompanying 
drawings in which: 

3 5 Figure 1 is a magnetic resonance spin echo image of 

the abdomen of a patient with liver iron 
overload acquired at an echo time of 6 ms 
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and a repetition time of 2500 ms; 



Figure 2 



shows a trace of maximal pixel signal 
intensities determined within a layer of 
subcutaneous fat surrounding the liver to a 
width of one pixel for the patient in 
Figure 1; 



Figure 3 



10 



shows the local maxima in signal intensity 
determined from the trace of the maximal 
pixel signal intensities in Figure 2, which 
correspond to the positions of the RF 
receiver points. 



15 Figure 4 



20 



illustrates the manner in which an angle Q n 
of a major axis of an elliptical contour 
to the decay in the RF signal intensity 
from a given RF receiver point is 
determined relative to the coronal plane 
from the trace of the maximal signal 
intensities throughout the subcutaneous 
fat; 



Figure 5A 



25 



shows a selected RF receiver point from 
Figure 3 and the two nearest neighbouring 
local minima; 



Figure 5B 



30 



35 



Figure 6 



shows the intensity of the local maximum at 
the RF receiver point and the intensity of 
the two nearest neighbouring local minima 
plotted as a function of distance in pixels 
from the RF receiver point and to which has 
been fitted an exponential signal intensity 
decay curve for the determination of R n ; 

shows the result of fitting the model 
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spatial intensity profile of the decay in 
signal intensity from the RF receiver 
points to the image intensities of selected 
pixels within the subcutaneous fat to 
5 obtain an estimate of the spatial variation 

in RF signal intensity throughout the liver 
of the patient in Figure 1; 

Figures 7A & 7B show the spin echo image of Figure 1 before 
10 (Figure 7A) and after (Figure 7B) reduction 

of the RF signal intensity non-uniformities 
throughout the liver via the estimated 
spatial variation in RF signal intensity of 
Figure 6 (with the contrast enhanced to 
15 accentuate image differences) ; 

Figure 8 shows a bi-exponential transverse 

magnetisation decay curve fitted to the 
signal intensity within the liver for a 

20 pixel p at position (x p ,y p ) of the patient 

in Figure 1 (for seven different spin echo 
image measurements) plotted as a function 
of the spin echo time, with the estimate of 
the initial signal intensity within the 

25 liver at zero echo time at pixel position p 

is included as an additional datum in the 
fit. Also shown is the expected signal 
decay rate within the subcutaneous fat 
relative to the RF spatial intensity 

30 profile value I ( (x p ,y p ) at pixel position 

P. 

Detailed Description of Preferred Embodiments 

35 As explained in greater detail below with reference to the 
accompanying drawings, in broad terms embodiments of the 
present invention provide a method of estimating the spatial 
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variation in RF signal intensity within magnetic resonance 
images of an object through the solution of a boundary value 
problem presented by a medium of essentially uniform spin 
density surrounding the object. From analysis of the 
5 surrounding medium, and the location of notional points of 
RF signal reception (hereinafter simply referred to as "RF 
receiver points") , a semi-empirical mathematical formulation 
of the decay profile of the RF signal intensity within the 
object is determined. The RF receiver points are located 

10 about the object in a manner consistent with the geometry of 
the RF detection coil elements in a magnetic resonance 
imager producing images of the object and the pattern of RF 
penetration throughout the object. These points are 

considered to be notional RF receiver points as they are not 

15 (usually) at the location of the actual RF detection coil 
elements of the imager but at locations which for the 
purposes of mathematical modelling have a substantially 
similar appearance to the positioning of the RF detection 
coil elements. The detection sensitivity of the RF receiver 

20 points is assumed to diminish in a concentric manner with 
increasing distance from the points, with the degree of 
ellipticity adjusted relative to the signal intensity 
pattern throughout the medium surrounding the object. To 
this end the RF receiver points may be located on the 

25 surrounding medium though this is not essential to the 
present invention. The sensitivity of each RF receiver 
point is also adjusted relative to the amplitude of the 
signal intensities throughout the bounding medium. The 
resulting semi-empirical mathematical formulation of the 

30 decay profile of the RF signal intensity throughout the 
object is then fitted to selected signal intensities from 
the medium bordering the object to obtain an estimate of the 
spatial variation in RF signal intensity within the object. 

35 More particularly, in the following description, a method is 
described for estimating magnetic resonance image 
intensities (as illustrated in Figure 6) within an object 
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such as a liver 10 (shown in Figures 1 and 7A) to map the 
spatial variation in RF signal intensity (Figure 6) within 
the object 10 that is due to instrumental and/or object 
effects. The method involves acquiring a magnetic resonance 
5 image of the object 10 surrounded by a medium of 
essentially uniform spin density, such as that provided by a 
layer of subcutaneous fat 12. A semi-empirical mathematical 
model is formulated of the spatial variation in RF signal 
intensity within the region encompassed by the surrounding 

10 medium 12 and which includes the object 10. Examples of two 
particular models are provided below, but in general terms, 
the model provides for a reduction, and in particular an 
exponential reduction, in RF signal intensity with 
increasing distance away from notional points of RF signal 

15 reception 14 (referred to as "RF receiver points 14") within 
the surrounding medium 12. The RF receiver points 14 are 
located at the local maxima of a selection of measured RF 
signal intensities from the surrounding medium 12. Prior to 
fitting, an estimate of the rate of decay in signal 

20 intensity from the RF receiver points is determined with 
reference to local minima 16 amongst the selected pixels 
from the surrounding medium 12. The model is then fitted to 
the image intensities of the selected pixels from the 
surrounding medium 12 to obtain an estimate of the spatial 

25 variation in RF signal intensity within the object 10 that 
is due to instrumental and/or object effects. 

Figure 1 illustrates an input image for the present method. 
In this regard Figure 1 illustrates a T2~weighted magnetic 

3 0 resonance image acquired of the abdomen of a patient in 
which the object under particular study is a liver 10. The 
layer of subcutaneous fat 12 which surrounds the liver 10 
provides a suitable surrounding medium of essentially 
uniform spin density which is used in the present method to 

35 estimate the spatial variation in RF signal intensity 
throughout the liver. In this instance, the image 

illustrated in Figure 1 shows a patient with liver iron 
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overload, the image being acquired at a spin echo time TE of 
6 ms and a repetition time TR of 2500 ms . 

Figure 2 illustrates a trace 18 of selected pixels within 
5 the subcutaneous fat 12 of the patient to be used in the 
estimation of the spatial variation in RF signal intensity 
within the liver 10. In one embodiment the selected pixels 
used to produce the trace 18 are pixels of maximal signal 
intensity 19 in respective radial line segments 17 of the 

10 image of the subcutaneous fat 12. Moreover the pixels 19 
are selected by an edge detection algorithm that locates the 
maximal signal intensities within the subcutaneous fat 12 to 
a width of one pixel around the periphery of the patient. A 
region of interest 20, drawn around the patient in Figure 2 

15 within which is shown the trace 18 of the maximal signal 
intensities within the subcutaneous fat 12, is used as a 
reference for the edge detection algorithm. The edge 
detection algorithm in essence look for the pixel of maximum 
intensity along each of the radial line segments 17 in the 

20 subcutaneous fat 12. It should be understood that the 
maximal intensities could be derived in relation to 
horizontal and vertical line segments from the region of 
interest 20 toward the centre of the body, rather than with 
reference to the radial line segments. 

25 

Figure 3 shows the local maxima 14 in signal intensity 
within the subcutaneous fat 12 determined from the trace 18 
of the maximal signal intensities shown in Figure 2. The 
local maxima 14 are located by sliding a window of constant 

3 0 width along the trace and for each position of the window 
ascertaining the maximum signal intensity for the length of 
the window. A second pass is then performed by sliding the 
window along the maximised signal trace and calculating the 
mean signal intensity of the maximised signal trace at each 

35 position of the window. Where the mean signal intensity is 
equal to the maximised signal intensity at a central window 
position, a local maximum 14 in signal intensity is 
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recorded. Ideally the window is selected to be at least 12 
pixels wide so that adjacent local maxima are no closer than 
12 pixels. 

5 As illustrated in Figure 3, seven local maxima 14i-14 7 are 
located. It would be recognised however that the number of 
local maxima will vary from patient to patient and object to 
object. The locations of the local maxima are taken as the 
positions of the notional RF receiver points for the 

10 mathematical description of the spatial intensity profile 
I(x,y) of the detected RF signal within the liver. While 
the local maxima in the bounding medium 12 are used as the 
RF receiver points 14, it is important to note that these 
points do not have to be located on the medium 12. The 

15 location of these points is used to determine the geometry 
of the spatial intensity profile. In providing the 

mathematical model of the spatial intensity profile I(x,y) 
it is assumed that signal strength diminishes or decays with 
distance from the points 14. Moreover, this decay is 

20 modelled as being an exponential decay and it is further 
assumed in the formulation of the model that RF signal 
intensity diminishes in a concentric manner from the points 
14. On the basis of these assumptions, the following two 
mathematical models of the spatial intensity profile I(x.,y) 

2 5 may be provided: 

Kx^tiy^'-"^ 3 "' 1 (i) 



Kx,y) = 2^in e (2) 

where x and y are image coordinates, N is the number of RF 
30 receiver points 14, I n is the signal intensity at the n th RF 
receiver point 14 n and is to be determined by the fitting 
procedure, R is the rate of decay in intensity of the 
detected RF signal with distance from the RF receiver points 
14 and is to be determined by the fitting procedure, x n is 
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the x position of the n RF receiver point 14 n , y n is the y 
position of the n th RF receiver point 14 n , f is the degree of 
ellipticity to the concentric reduction in RF signal 
intensity with increasing distance from the RF receiver 
5 points 14 and is to be determined by the fitting procedure, 
and Q n is the angle of the major axis of the elliptical 
contour to the decay in the RF signal intensity from the n th 
RF receiver point 14 n with the coronal plane. In relation to 
Figure 3, the intensity values of the local maxima are taken 
10 as initial estimates of the signal intensities i n at the N 
RF receiver points throughout the subcutaneous fat. For the 
example under discussion, N = 7. 

The mathematical model (2) above is a more general form of 
15 mathematical model (1) . Both models provide for an 

exponential decay in signal intensity or strength from the 

RF receiver points in an elliptical manner. However model 
(2) also allows for orientation of the major and minor axes 

of the elliptical decay in signal intensity for each of the 
20 RF receiver points 14. The model (2) reduces to the model 

(1) if it is assumed that major axis of the elliptical decay 

contours is aligned with the coronal plane 13, this is 

equivalent to Q n being 0°. 

25 Figure 4 illustrates the manner in which the angular 
alignment of the major axis of the elliptical contour of the 
decay in the RF signal intensity from a given RF receiver 
point 14 is determined relative to the coronal plane 13 from 
the trace of the maximal signal intensities throughout the 

3 0 subcutaneous fat as shown in Figure 2. A window kernel 15 
is centered on a given RF receiver point and the two most 
extreme positions at which the trace 18 of the maximal 
signal intensities passes through the perimeter of the 
window kernel 21 are determined. For the example in Figure 

35 4, the window kernel 15 is square with a width of 25 pixels. 
The angle Q of the line bisecting the two said extreme 
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positions on the perimeter of the window kernel is then 
determined relative to the coronal plane 13, as denoted in 
Figure 4. For the n th RF receiver point 14 n , this angle is 
denoted Q nr and is used in the model (2) of the spatial 
5 intensity profile I(x,y) of the detected RF signal intensity 
within the liver. 

Figures 5A and 5B illustrate how an initial estimate of the 
rate R of decay in intensity of the RF signal with distance 

10 from the RF receiver points 14 is determined by way of 
example to one RF receiver point 14 x . Figure 5A shows a 
selected RF receiver point 14i from Figure 3, and the two 
nearest neighbouring local minima 16a and 16b, determined in 
a manner consistent with that for location of the local 

15 maxima in the subcutaneous fat 12. Figure 5B shows the 
intensity of the local maxima at the RF receiver point 14 
and the intensity of the two neighbouring local minima 
plotted as a function of distance in pixels from the RF 
receiver point 14. An exponential decay curve is fitted to 

20 the signal intensity data points with the rate R n of decay 
in signal intensity as the unknown parameter, and with the 
intensity at zero distance from the RF receiver point 14i 
fixed to the intensity of the local maxima at that point. 
This procedure is repeated for each RF receiver point 14, 

25 for which the estimates of the rate R n of decay in signal 
intensity with distance from each RF receiver point 14 are 
averaged to give an initial estimate of R, the rate of decay 
in intensity of the detected RF signal with distance from 
the RF receiver points. 

30 

Figure 6 shows the result of fitting the model spatial 
intensity profile I(x,y) to the selected image intensities 
from the subcutaneous fat 12 to obtain an estimate of the 
spatial variation in RF signal intensity throughout the 
35 liver 10. The accuracy of the fit is improved by removing 
any signal intensities from the trace 18 of the maximal 
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signal intensities throughout the subcutaneous fat 12 that 
are either inconsistent with the presence of fat or are 
otherwise erroneous. The projected signal intensities may be 
used to rescale the magnetic resonance image intensities 
5 within the liver and thus minimise the spatial variation in 
RF signal intensity throughout the liver, as shown in 
Figures 7A and 7B. Figure 7A is the spin echo image of 
Figure 1 redisplayed with increased image contrast, and for 
which a region of interest has been drawn about the liver 

10 10. Figure 7B shows the liver signal intensities with a 
reduced variation in RF signal intensity, in which the liver 
signal intensities have been divided through by a RF spatial 
attenuation profile determined from the fitted spatial 
intensity profile I(x,y). In this instance, the RF spatial 

15 attenuation profile was determined by dividing the fitted RF 
signal intensity profile I(x,y) by the minimum value of 
I(x,y) within the region occupied by the liver. 

Figure 8 illustrates application of the estimated spatial 

20 variation in RF signal intensity within the liver to 
facilitate the calculation of transverse relaxation rate 
(R 2 ) within the liver. This application is particularly 
useful where the transverse relaxation rates within the 
liver are determined from a series of spin echo images where 

25 the signal intensities are predominantly less than 50% of 
the initial signal intensity at zero echo time. The RF 
signal intensity profile I(x,y) fitted to the maximal signal 
intensities within the subcutaneous fat is transformed into 
an estimate of the image intensities within the liver at 

30 zero echo time by dividing I(x,y) by the percentage 
remaining signal intensity expected within the subcutaneous 
fat at the given echo time, and then dividing this result by 
the ratio of the hydrogen proton spin density expected of 
the subcutaneous fat relative to the hydrogen proton spin 

35 density expected of the liver. The estimate of the image 
intensities within the liver at zero echo time may be then 
included as additional data points at zero echo time in 
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calculation of the transverse relaxation rates within the 
liver. This is illustrated in Figure 8 where the signal 
intensity within the liver for a pixel p at position (x p ,y p ) 
from seven separate spin echo time measurements is plotted 
5 versus echo time, with the estimate of the initial signal 
intensity within the liver at zero echo time (determined 
from the 6ms spin echo image) at pixel postion p is included 
as an additional datum at zero echo time. In determining the 
initial signal intensity within the liver, the percentage 

10 remaining signal intensity expected within the subcutaneous 
fat at the given echo time of 6ms as shown in Figure 8 is 
effectively given by the values X divided by Y which are 
marked on the signal intensity axis. The value for the 
percentage remaining signal intensity expected within the 

15 subcutaneous fat is in practice estimated from the expected 
signal decay rate for the subcutaneous fat (R2Fat) relative 
to the given echo time (TE) by the expression exp (-R2 F at . TE) , 
where exp ( ) is the exponential function. The ratio of the 
hydrogen proton spin density expected of the subcutaneous 

20 fat relative to the hydrogen proton spin density expected of 
the liver is evidenced in Figure 8 by the signal intensity 
axis values of Y divided by Z. The transverse relaxation 
rate characteristic of the liver in the example of Figure 8 
was determined by fitting a bi-exponential signal decay to 

25 the data points. The inclusion of the estimated signal 
intensity at zero echo time implicitly improves the accuracy 
of the calculated R2 by including extractable information 
about signal intensities at points in time much earlier than 
that captured in the shortest echo time image. This method 

30 is especially important in the accurate modelling of bi- 
exponential transverse magnetisation decay processes. 

From the above description, it will be apparent to those 
skilled in the relevant arts that the present method has 
35 numerous advantages and benefits over prior art methods for 
correcting spatial non-uniformities in RF signal intensity 
within magnetic resonance images. Specifically, the present 
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method enables the determination of a RF signal intensity 
correction surface that is particularly sensitive to non- 
uniformities in the RF receiving coil profile as well as 
signal intensity attenuation within the object under study. 
5 This is facilitated through the concurrent imaging of a 
medium of essentially uniform spin density bordering the 
object from which the pattern of RF signal intensity 
attenuation within the object is estimated. A semi-empirical 
mathematical model of the spatial variation in RF signal 

10 intensity within the object is determined through the 
location of notional RF receiver points about the boundary 
medium, from which a concentric decay to the detected RF 
signal intensity is assumed with increasing distance from 
the said points. The concentric pattern of decay from each 

15 RF receiver points is taken to be elliptical, with the major 
axes of the elliptical contours approximately aligned with 
the RF coil. The ellipse is a reasonable model of the decay 
profile in RF signal intensity that results from the 
convolution of the RF coil geometry with the internal 

20 attenuation effects from the object under study. 

Now that an embodiment of the invention has been described 
it will be apparent to those skilled in the relevant arts 
that numerous modifications and variations may be made 

25 without departing from the basic inventive concepts. For 
example, the RF receiver points need not be constrained to 
location at the maximal signal intensities on the boundary 
medium, or even on the boundary medium itself. Further, the 
decay in signal intensity with increasing distance from the 

30 RF receiver points need not follow an exponential decay, or 
the same decay rate from each point. In addition, the signal 
intensities within the object under study may be rescaled by 
direct division with the fitted signal intensity profile to 
form a signal intensity ratio of the object image 

35 intensities relative to the boundary medium intensities. For 
example, in abdominal imaging, the signal intensity ratio of 
liver to subcutaneous fat may formulated at every point 
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within the liver with reference to the projected signal 
intensities determined from the subcutaneous fat. Finally, 
it * should be noted that the method described in this 
specification is broadly applicable to the estimation of 
5 magnetic resonance image intensities in which there is a 
continuous decay in signal intensity with distance from 
several points of maximal signal intensity within the image. 

All such modifications and variations that would be obvious 
10 to a person of ordinary skill in the art are deemed to be 
within the scope of the present invention the nature of 
which is to be determined from the above description and the 
appended claims. 
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